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Abstract 

We investigate perturbed monomial dynamical system over F p given 
by iterations of i Hi" + c modp, where c g F p . Instead of study the 



systems one at a time we study all of them at the same time. The complex 



distibution of periodic points is visualized in the so called Periodic Point 
Diagram, which can be seen as a discrete version of the classical Bifurca- 
tion Diagram. We also prove some general results about the distribution 
of periodic points. We end the article with a conjecture about the total 
number of periodic points. 

1 Introduction 

^-i- ' We consider the dynamical system given by iterations of 

•■ h c (x)=x n + c (1.1) 

_~ , over the prime field ¥ p , the set {0, 1, . . . , p — 1} with addition and multiplication 

modulo p. We let n > 2 be an integer and c <E {0, 1, . . .,p — 1}. By h r c we 
mean the r-fold composition of h c . The dynamics of h c (x) changes dramatically 
with the value of the parameter c. Classically we have a bifurcation if there is 

j^ , a sudden change in the dynamics at a certain parameter value, see for example 

[5]- If we use the classical definition of a bifurcation point we can say that we 
have a bifurcation for every value of c. When describing the dynamics we will 
concentrate on the number of periodic points. A point a is said to be a periodic 
point if h r c (a) = a for some positive r, the smallest such r is called the period 
of a. We say that a is an r-periodic point. Since we are dealing with a finite set 
(all residues modulo a prime p) every point will after a number of iterations be 
mapped onto a periodic point. The set of points can therefore be partitioned 
into two sets, the periodic points and the preperiodic points that will eventually 
be mapped onto periodic points. 

We will use the notation F* for the multiplicative group of ¥ p . We will 
visualize the dynamics of h c on ¥ p by a directed graph Graph(V, E), where the 
the vertex set is V = ¥ p and the edge set E, is the set {(x, h c (x)),x £ V}. 



Example 1.1. Let us study the dynamics of the h c (x) = x 2 + c on F7 for all 
seven possible choices of c. We iterate the map x »->• x 2 + c mod 7. The result 
is shown in Figure 11.11 Here we can see that the dynamics changes a lot if we 
change the value of c. 

The systems x M- x 2 + c on ¥ p have been studied before. In [8] Rogers 
discribes the connected components of the directed graph of a; i-> x 2 . Rogers 
stated that for c ^ the dynamics seems to be beyond description. In [11] Vasiga 
and Shallit succeeded in describing the dynamics of x i-> x 2 — 2. Also Gilbert et 
al in [2] contributed to the understanding of the dynamics of x 1— *• x 2 — 2. We 
can see this already in Figure 11.11 the graphs for c = 5 and c = seems more 
symmetric than the other graphs. In [3] the number of cycles (periodic orbits) of 
x »->• x n mod p are calculated by using Mobius inversion formula. Asymptotical 
behavior for the number of periodic points when p — >■ 00 where investigated. 
The results from [4] where extended in [6] and in [7]. 

In this article we will construct a bifurcation diagram analogously to the 
classical diagram describing for example the bifurcations of the logistic map 
x M> rx(l — x) over the real numbers. See for example [3J for details on this 
dynamical systems. In Figure 11.21 a bifurcation diagram for the logistic map 
is given. A variant of the Mathematica program suggested by J3J was used for 
producing the picture. 

This article is organized as follows: In Section [5] we introduce the Periodic 
Point Diagram (PPD) and compare it with the classical bifurcation diagram. 
We also prove some general properties of the PPD. Then in Section [3J and 
Section @] we investigate the properties of PPD in more detail for n = 2 and 
n > 2, respectively. In Section [5] we formulate a conjecture on the asymptotical 
behavior of the total number of periodic points of set {x M> x" + c; c € ¥ p } of 
dynamical systems when p — >• 00 . 



2 The periodic point diagram 

First of all there is an important difference between the bifurcation diagram 
we will construct for x n> x n + c mod p and the classical one; our diagram is 
discrete. We will include every possible value of both x and c. We construct the 
diagram in the following way: On the horizontal axis we have all the possible 
values of c, c G F p and on the vertical axis we have the possible values of x, 
x <e¥ p . For each c we mark the x- values that are periodic points. 

Example 2.1. Let us consider the dynamical system x 1— > x 2 + c mod 7. In 
Figure 12.11 the bifurcation diagram for this system is shown. It is hard to see 
any pattern from one c to the next. However, there seems to be diagonal lines 
with no periodic points. 

The diagram in Figure [2~T1 was made by using a Mathematica program, that 
naively computes the periodic points for each c: Iterate every point in ¥ p , p 
times then start saving the next p iterates. Then take the union of all resulting 
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Figure 1.1: The dynamics of the seven systems x i— > x 2 + c mod 7. 
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Figure 1.2: The bifurcation diagram of x M> rx(l — a;) for r s [2, 4]. This is just 
iterations of the starting value x = 0.5. 
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Figure 2.1: The bifurcation diagram of x n- x 1 + c mod 7. 
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Figure 2.2: The diagonal line j = i + 3 mod 7 in the PPD(/i c , 7). 



points. In this way we are sure to have found all the periodic points and no 
preperiodic points. Of course there a more effective algorithms that uses for 
example Floyd's cycle finding algorithm (Tortoise and Hare), see p] [5], but 
creating the diagrams in this article the naive way is enough. 

We now make some formal definitions. First of all, instead of saying "bifur- 
cation diagram" we call the diagram Periodic point diagram (PPD), it seems 
like a better name since we do not really have the same kind of bifurcation in 
our finite case as in the classical case. The PPD of the parametrized dynamical 
system given by iterations of the mapping x i— > h c (x) modulo p is denoted by 
PPD(h c ,p). Let a G ¥ p . By a diagonal line in the PPD(h c ,p) we mean all 
points (i,j) in the diagram satisfying the equation j = i + a mod p. For a = 
we have the diagonal from the lower left corner to the upper right corner. All 
other diagonal lines are divided into two parts, see Figure 12.21 



Definition 2.2. By a desert line of PPD(h c ,p) we mean a diagonal line with 
no periodic points. 

Example 2.3. Let us now consider the dynamical system x t-» x 2 + c mod 71. 
In Figure |2~31 we have its PPD. Here the desert lines are clearly visible. 



Before we enter the discussion on how many desert lines there are, we take 
a look at some conditions for h c to be bijective. 
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Figure 2.3: The bifurcation diagram of x i— > x 2 + c mod 71. 

Lemma 2.4. // /io is bijective then h c is bijective for all c € {0, 1, . . . ,p — 1}. 

Proof. Assume that the map ho(x) — x n is bijective. Let x, y G ¥ p . li x n + c = 
y n + c then x" = y n and this implies x = y since x n- a;™ mod p is bijective. D 

In fact, we can generalize the lemma to: If h c is bijective for one c it is 
bijective for all. Moreover, if the map h c is bijective then the dynamical system 
can't have any preperiodic points. Hence the dynamics has a pure cycle structure 
if and only if h c is bijective. 

Theorem 2.5. The dynamics of h c {x) — x n + c has a pure cycle structure if 
and only if gcd(n,p — 1) = 1. 

Proof. In this proof we will use results from elementary number theory. See for 
example [9] for details on primitive roots and solutions of congruence equations. 
The theorem follows if we can prove that i4i" mod p is bijective if and only 
if gcd(n,p — 1) = 1. Let p be a primitive root modulo p. Let x,y 6F* then 
there are i, j € {0, 1, ... ,p — 2} such that x = p? and y — p? . We have x n = y n 
if and only if ni = nj (mod p — 1). This equation has the unique solution i = j 
(mod p — 1) (hence i = j since they are both in {0, 1, . . . , p — 2}) if and only if 
gcd(n,p — 1) = 1. Also note that x n = in F p if and only if x = 0. □ 

In the next section we will prove that there are exacly (p — l)/2 desert lines 
in PPD(.T 2 + c,p) forp> 3. 



3 The quadratic case 

In this section we look at the PPD for the quadratic mapping x t-t x 2 + c mod p. 
We first recall the definition of quadratic residues and quadratic non-residues. 

Definition 3.1. Let p be a prime and let a be an integer such that gcd(a,p) = 1. 
We say that a is a quadratic residue if the equation x 2 = a has a solution in ¥ p . 
If it doesn't have a solution, then a is called a quadratic non-residue. 

We have the following well known result for quadratic residues, see [3J: 

Theorem 3.2. Let p be an odd prime. Then there are (p — l)/2 quadratic 
residues and (p — l)/2 quadratic non-residues. 

Theorem 3.3. Let p > 3 be a prime. If a is a quadratic non-residue modulo 
p then j = i + a mod p is a desert line. If a is a quadratic residue there are 
two fixed points on the line j = i + a mod p. If a = there are one fixed point 
on the line j = i (mod p). Hence, there are exactly {p — l)/2 desert lines in 
PPD(ft c ,p). 

Proof. We have an r-periodic point at position (i,j) in the PPD-diagram if 
h\(j) — j mod p. Consider the points (i, j) on the diagonal line j = i + a mod p. 
Assume that we have an r-periodic point on this line. That is, for some j £ 
{0, 1, ... ,p — 1} we have 

h\(i + a) = i + a 

We have 

(hi' 1 (i + a)) 2 + i = i + a, 

and 

(hr\i + a)f=a. (3.1) 



Hence, if a is a quadratic non-residue then (|3.1[) has no solution and j '• — i + 
a mod p is a desert line. From Theorem 13.21 it follows that we have at least 
(p — l)/2 desert lines in the PPD. 

Let r = 1 and let a be a quadratic residue. We have from (|3.1[) that (i+a) 2 = 
a. There are exactly two values of i that solves this equation since there are 
exactly two square roots of a quadratic residue. For a — there is exacty one 
fixed point (for i = 0). So, we have exactly (p — l)/2 desert lines. □ 

4 The n-power case 

In this section we will generalize our investigations from last section to the PPD 
of x <— > x n + c mod p. First we note that if n and p — 1 are relatively prime then 
it follows from Theorem 12.51 that we have a pure cycles structure for all values 
of c. This means that we have no desert lines in this case. Moreover, the whole 
PPD is filled, since every point is a periodic point. Before we learn more about 
the desert lines we recall some facts about n-power residues. 



Definition 4.1. Let p be a prime and let a £ F*. Then a is said to be an n- 
power residue if x n — a has a solution in ¥ p . Otherwise, a is called an n-power 
non-residue. 

Theorem 4.2. We have that a S ¥ p is an n-power residue if and only if 

a (p-!)/ gcd(n.p-l) _ ^ 

There are (p — l)/gcd(n,p — 1) such a in ¥ p . For each n-power residue a the 
equation x n = a has gcd(n,p — 1) solutions in ¥ p . 

See for example |10) for a proof of this theorem and for more details on 
n-power residues. 

Theorem 4.3. We have that j = i + a is a desert line in PPD(h c ,p) if and 
only if a is an n-power non-residues modulo p. Hence there are exactly p — 1 — 
(p — l)/gcd(n,p — 1) desert lines. Moreover, if a is an n-power residue then 
there are gcd(p — 1, n) fixed points on the line j — i + a. For a — there is one 
fixed point (for i = 0). 

Proof. Consider the points (i,j) on the diagonal line j = i + a. Assume that we 
have an r-periodic point on this line. That is, for some i G {0, 1, . . .,p — 1} we 
have 

h\{i + a) = i + a. 

We have 

(h\~ 1 (i + a)) n +i = i + a, 

and 

(hr\i + a)) n =a. (4.1) 

Hence, if a is an n-power non-residue then j = i + a is a desert line. 

If r = 1 and a is an n-power residue then from (|4.1[) we have (i+a)"+i = i+a 
(mod p) and hence (i + a) n = a. From Theorem 14.21 it follows that there are 
gcd(n,p — 1) different i £ FI satisfying this equation. For r — 1 and a = 0wc 
have the unique solution i — 0. Hence there are exactly 

p — 1 p — 1 

V- 1 ^7 7T = —77 rr(gcd(n,p - 1) - 1) 

gcd(n,p-l) gcd(n,p-l) 

desert lines in the PPD. □ 

Example 4.4. Let us look at the PPD of the dynamical system x n> a; 5 +c mod 
71 in Figure I4T41 Here we see the increased number of desert lines compared to 
the x n- x 2 + c mod 71 in Example 12.31 



We end this section with an observation regarding symmetry when n is odd. 

Definition 4.5. By the reduced Power Point Diagram of h c (x) over F p we mean 
the diagram we get by removing the column c = and the row x — from the 
PPD. 
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Figure 4.1: The PPD(x 5 + c, 71) 



Theorem 4.6. If n is odd then the reduced PPD of h c over F p is symmetric 
with respect to the line j = —i mod p. That is, if a is an r-periodic point of 
h c (x) then —a is an r-periodic points of h- c (x) (—a and —c are the additive 
inverses of a and c modulo p). 

Proof. We first observe that h- c (— a) = (—a) n —c = —h c (a) and that h r _ c (~a) = 
—h r c (a) for all r > 1. Assume now that a is an r-periodic point of h c (x). Then 

h r _ c (-a) = -h r c {a) = -a. 

Moreover, h_ c (—a) ^ —a for all d < r since hf(a) ^ a for such d. Hence —a is 
an r-periodic point of h- c {x). □ 

5 Conjectures about the number of periodic points 



From Theorem 14.31 it follows that there are exactly p fixed points in total to the 
set \h c {x) : c £ F p } of perturbed monomial systems modulo p. Hence, among 
the marked points in the PPD there are exactly p that are fixed points. 

In this section we will consider the problem of computing the number of 
periodic points for all values of c, that is the number of marked points in the 
PPD. Let Per(/i c ,p) be the number of periodic points of the class of dynamical 
systems x i— > h c (x) mod p (the number of marked points in the PPD(ft, c ,p)). 
The only case when it is easy to compute Per(/i c ,p) is when gcd(rt,p — 1) = 1. 
Then we have Per(/i c ,p) = p 2 . From l4~3l we have the following rough estimate 

p<Per(h c ,p)< ^~ 1} 1V (5.1) 

gcd(ra,p-l) 

We have used the fact that there can be at most p periodic points on each 
non-desert line. 

We will use a Mathematica program for computing the number of periodic 
points. The program is based on Floyd's cycle finding algorithm [I] [5]. The 
naive program described in Section [5] for producing the PPD:s is too slow for 
large primes. Let us look at some examples. 

Example 5.1. Let us calculate the number of periodic points of h c (x) = x 2 +c 
(mod p) for the first 1000 primes. The result is shown in Figure [5TT1 

Example 5.2. Consider the dynamical system h c (x) = x 12 + c. In Figure 15721 
Per(x 12 + c,p) is plotted for the first 1000 primes. Here we see that the diagram 
seems to contain four different curves. These are branches corresponding to the 
four possible different values of gcd(12,p — 1) for odd primes p. The possible 
values are 2, 4, 6 and 12. The lowest branch corresponds to 12 and the highest 
to 2. 

We call the kind of diagrams shown in the examples Total Periodic Diagrams 
(TPD). The Total Periodic Diagram for h c (x) = x n + c for the first m primes 
is denoted by TPD(ft, c ,m). 
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Figure 5.1: Per(x 2 + c,p) for the first 1000 primes. Note that the value of the 
horizontal axis is the ordinal of the prime number. 
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Figure 5.2: Per(cc 12 + c,p) for the first 1000 primes. 
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Conjecture 5.3. Let n > 2 be an integer. Let h c (x) = x n + c mod p. The 
TPD(/i c ,to) then contains as many branches as there are possible values of 
gcd(n,p- 1). 

It would of course also be interesting to find a simple asymptotic estimate 
for Per(a; n + c,p) when p — > oo. Thus, we are interested in finding a function 
fd(p) such that 

Per(x n + c,p) 
hm —— = 1, (5.2) 

gcd(n,p— 1)— a 

for a possible value d of gcd(n,p — 1). 

6 Discussion 

Some of the results of this article can be generalized to arbitrary functions f(x) 
defined on F p . Desert lines will occur also in the general case. There will be as 
many desert lines as there are a £ ¥ p such that f(x) — a has no solutions in F p . 
It is also possible to extend some of the results to more general structures than 
the fields F p , for example the rings of integers modulo an composite integer. 
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